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Highly focused supersonic microjets: 
numerical simulations 
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By focusing a laser pulse inside a capillary partially filled with liquid, a vapour bubble 
is created which emits a pressure wave. This pressure wave travels through the liquid 
and creates a fast, focused axisymmetric microjet when it is reflected at the meniscus. 
We numerically investigate the formation of this microjet using axisymmetric boundary- 
integral simulations, where we model the pressure wave as a pressure pulse applied on 
the bubble. We find a good agreement between the simulations and experimental results 
in terms of the time evolution of the jet and on all parameters that can be compared 
directly. We present a simple analytical model that accurately predicts the velocity of 
the jet after the pressure pulse and its maximum velocity. 
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1. Introduction 



In recent experiments by Tagawa et al. (2012), it was found that microscopic jets that 



travel at a speeds up to 850 m/s can be created by focusing a laser pulse inside a liquid- 
filled capillary that is open at one end. Besides the high velocity, the jets were found to be 
highly reproducible and controllable. The laser pulse used in the experiments, which has 
an energy of the order of 100 /zJ, results in the formation of a vapour bubble accompanied 



by a pressure wave ( Bell & Landt 1967 Felix & Ellis 1971 ). This pressure wave is reflected 



at the free surface, where the jet is formed. The shape of the free surface was found to 
play a crucial role in the formation of the jet, as it is responsible for focusing the liquid 
into a jet. 

In this paper, we present numerical simulations which accurately reproduce the evo- 
lution of the shape and the velocity of the jets observed in the experiments described 
in Tagawa et al. (2012 1. We use axisymmetric boundary integral simulations where we 



model the effect of the traveling pressure wave by applying a short pressure pulse on a 
bubble with a constant amplitude Ap and a (short) duration At such that the resulting 
impulse per unit area ApAt is of the order of 10 Pa • s. Figure [T] shows the formation of a 
jet in the experiment together with a result from our boundary integral simulations. We 
provide a theoretical analysis that gives the correct scaling for the jet speed as a function 
of the contact angle, bubble distance and pressure pulse. 

Jets that are formed from a meniscus have been studied in several cases with different 



sources for the driving pressure: In Duchemin et al. (2002), the pressure was provided 
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Figure 1. Jet formation in experiment (background images) and simulation (black and white 
solid lines). Absorbed laser energy was 365 fiJ, distance between laser focus and meniscus was 
600 (jm, the tube radius is 250 /jm. Pressure amplitude used in the simulation was 1581 bar, 
pressure duration was 50 ns, and the initial bubble radius was 25 /jm. 



by surface tension, and Antkowiak et al. (20071 created a pressure pulse by impacting a 



test tube on a hard surface. In Lorenceau et al. ( 2002 ) and Bergmann et al. ( 2008 ) the 



driving was provided by hydrostatic pressure. In all these cases the scale of the jets is 
much larger and the velocities are at least an order of magnitude smaller than those that 
we study here. 

The formation of the jet is different from the jet that follows from the collapse of a 
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liquid void (Longuet-Higgins & Oguz 1995 Hogrefe et al. 1998 Bergmann et al. 2006 



Gekle et al. 



2009), where the jet initiates from a geometric singularity. In that case the 



size and the initial speed of the jet following the collapse of a cavity is therefore not set 
by the typical size and velocity of the experiment. Neither does the theory of a hyperbolic 
jet (Longuet-Higgins 1983) apply here. The main ingredient of the latter is a hyperbolic 
radial inflow from infinity, which in our setup is impossible due to the confinement of the 
tube. In fact, the jet that we study scales with the size of the capillary, and the speed is 
controlled by a combination of driving and geometry. 

The paper is organized as follows. First we introduce the numerical method in §[2j 
including a discussion of the geometrical setup and a discussion of how the simulations 



connect to the experiments of |Tagawa et al. (20121). Subsequently, we present the results 
from the numerical simulations and compare these to the experiments in §[3j discussing 
the influence of the various parameters on the maximum velocity of the jet. In §|4]wc 
then derive an analytical model that predicts the achieved velocities, and we end with 
the conclusions and discussion (§[5]). 



2. Numerical setup 

We perform numerical simulations using a boundary integral code, as described in Oguz 
fc Prosperetti] (|1993[ ), |Bergmann et al.\ ( |2009[ ), |Gekle fc Gordillo| ( |2010a[ ) and |Gekle k 
Gordillo (20106 ). Here, we repeat the basic principles and methods, and elaborate on the 



parts that are specific for our case. 

We approximate the flow in our system to be incompressible, inviscid and irrotational, 
so that we can describe the velocity field v as the gradient of a potential 



which satisfies the Laplace equation 



v = V</» 



= 0. 



(2.1) 



(2.2) 



Using Green's identity, the potential at any point inside the liquid domain can be de- 
scribed by an integral over the boundary containing (j) and <j> n , where (j) n is the spatial 
derivative of <j> in the direction normal to the boundary. The system can be solved if 
at every point on the boundary either <j) or <p n is known. Solving the system is greatly 
simplified by imposing axial symmetry, reducing the surface integrals to line integrals. 
This simplification is justified by the axial symmetry observed in the experiments. On 
stationary solid boundaries we have <p n = 0, and on the free surface we know the potential 
after time-integrating the unsteady Bernoulli equation 



Of 
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2 p 



gz 



(2.3) 



with Ap — p g —p a the pressure of the ambient vapour minus the atmospheric pressure, k 
the curvature, a the surface tension, p the liquid density, g the gravitational acceleration 
and z the vertical coordinate. Due to the size and the time scale of the experiment, 
the gravitational component gz can be neglected. After solving the boundary integral 
equation, <f> and <p n are known everywhere on the boundary, and the new position of the 
free surface can be achieved by time-integrating the kinematic boundary condition 

| = V * (2.4) 

Due to the absence of viscosity, some form of surface smoothing is necessary to keep the 
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Figure 2. The axisymmetric numerical setup. A tube of radius Rt is filled with liquid. The 
liquid-air interface has the shape of a spherical cap, with a contact angle 0. We position a 
bubble with a radius of 1/10 of the tube radius at a distance A from the meniscus. The dashed 
line in the center represents the axis of symmetry. 



simulation stable. We use the node-shifting technique described by |Oguz &, Prosperetti| 



(1990), according to which new nodes are placed half way between all existing nodes, 
after which the original nodes are removed. This method effectively removes instabilities 
that are related to the node spacing everywhere on the free surface, except at the node 
on the symmetry axis because this node cannot be shifted or removed. We found that 
in our situation this node was subject to these instabilities, and therefore we applied an 
additional smoothing to it. This was done with the help of quadratic extrapolation of 
the position and the potential, using the two nodes next to the axis of symmetry and the 
symmetry condition. We verified that the numerical solutions were not sensitive to the 
amount of smoothing that we applied. 

Because we are investigating a liquid inside a capillary, we have to take into account 
a moving contact line. The node that connects the liquid surface to the solid boundary 
can be considered both as part of the free surface and of the capillary wall. In solving the 
boundary integral equation, we treat this connecting node as a the latter, and impose 
<fi„ = on it. Implementation of the actual moving contact line with a dynamic contact 



angle as described by Voinov ( 1976 ) would not be appropriate here since this is based on a 
balance between surface tension and viscosity, while our simulations are inviscid. Instead, 
we calculate the new position of the connecting node by extrapolating the nodes next 
to the connecting node. This method is similar to the one used by|Oguz &: Prosperetti| 



(1993), where the node connecting to the solid was displaced so that the contact angle 
remained at 90° . The exact implementation turns out only to have a non- negligible effect 
only close to the contact point and not to be important for the development of the jet. 
Comparing the extrapolation method and fixed contact angles between 60° and 120° 
resulted in less than 1% variation of the maximum jet velocity. 

2.1. Initial condition 

The computational domain that we use is closed at one side in the shape of a half-sphere; 
the free surface is at the opposite side of the tube (see figure [2J). Because the diameter 
of the tube is much smaller than the capillary length, gravity can be safely neglected, 
and the free surface adopts the shape of a spherical cap. The initial shape of the free 
surface in the simulations can therefore be defined using only the static contact angle 9. 
Note that the contact angle 9 only serves as an initial condition, and we do not impose 
a dynamic contact angle. A bubble with a radius 1/10 of the tube radius is positioned 
at a distance A from the free surface. The distance of the bubble from the closed end of 
the tube does not have an influence on the simulations, as long as it is a few tube radii 
or more. 
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2.2. Pressure wave model 

In the experiment, a pressure wave is created by vaporizing a small amount of liquid 
with a laser pulse. This abrupt vaporization is responsible for a very large increase in the 
pressure in a small volume, which results in a pressure wave that travels through the tube 
and reflects on the free surface (Tagawa et a/.|20"l2). As argued in that paper, the initial 



velocity Vq of the free surface is connected to the pressure wave strength Ap sa ^pcVo, 
where c is the speed of sound in the liquid. E.g., for Vq ~ 10 m/s and c = 1497 m/s wc 
find Ap ~ 75 bar. 

The reflections of the pressure wave on the free surface and the wall of the tube 
ultimately result in a pressure gradient between the vapor bubble and the free surface, 
so that the entire liquid volume in between will start to move. Starting from this very 
early point in time, the dynamics of the system is expected to be well described by the 
potential flow boundary integral model employed in this paper. In the simulation, we 



model this pressure wave by applying a pressure pulse on the bubble ( Ory et al. 2000 ) . 
This pressure pulse has a typical amplitude Ap and a duration At. Figure [3j shows that 
if At is small enough, the only relevant value is the product ApAt (which is in the order 
of 10 Pas) and the pressure pulse can be assumed to be instantaneous. For simplicity, 
we will keep At at 50 ns, and vary only the pressure amplitude. With this choice, the 
pressure amplitude easily reaches values in the order of 10 3 bar, i.e., in excess of the 
critical pressure of water (« 220 bar). This does not need to worry us too much, since 
not Ap itself but the product ApAt determines the course of events. After the pressure 
pulse, the pressure inside the bubble is set to to account for the rapid condensation of 



the vapour in the bubble as was also done by Ory et al. (2000). A more sophisticated 



model where the heat exchange is taken into account for the growth and collapse of a 



vapour bubble created by a laser pulse can be found in Sun et al. (2009). In our case, 
applying a perfect gas law as well as incorporating heat transfer only resulted in marginal 
differences in the jet. We therefore use a simpler model here, which minimizes the number 
of unknown adjustable parameters. The good agreement between the simulations and the 
experiments convinces us that our model is accurate enough to describe the physics that 
create the jets observed in experiments. 



3. Numerical results 

We will show here the results from the numerical simulations, compare them to the 



experiments of Tagawa et al. (2012), and discuss the influence of various parameters. 
More specifically, we will first compare the time-evolution of the shape and speed of the 
jet from the numerical simulations with the experiments. After this, we will focus on the 
maximum velocity of the jet by investigating the influence of the contact angle, surface 
tension, pressure impulse, tube radius, and bubble distance, and discuss these results in 
the context of the experiments. 

In figure [T] we find 11 snapshots from the experiment, overlaid with the corresponding 
BI results. At t = the laser pulse is absorbed which is modelled in the simulation 
as a bubble pressurized to 1581 bar for a time span of 50 ns. The initially concave 
interface starts to move instantly and goes through an almost flat stage at t = 7 /is to 
the development of a jet (t = 14 — 56 /is). The maximum jet speed is reached between the 
second and the third frame at t = 10 /xs, when the jet just starts to form. The position of 
the tip of the jet is fully reproduced, as well as the overall shape of the jet. The bubble 
in the experiment is created near the wall of the capillary, resulting in a difference in 
shape, but its size is reproduced by the simulations at least up to t = 35 /is, as can 
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Figure 3. The velocity of the tip of the jet as a function of time for different values of At. The 
product AtAp is kept constant at 15.20 Pa ■ s, showing that the velocity of the jet only depends 
on this product when At is small enough. The results for At = 25 ns (circles, Ap = 6080 bar) 
and At = 50 ns (solid line, Ap = 3040 bar) overlap, whereas a significantly larger duration of 
the pressure pulse At = 6.4 fjs (dashed line, Ap = 23.75 bar) results in a different velocity and 
a different evolution of the velocity. 



be appreciated by considering volume conservation in the system. Because the full free 
surface of the experimental jet is reproduced by the simulation up to t = 35 /is, we know 
that the volume of the bubble in the simulation is also the same as in the experiment. 
After this there is a slight difference in the collapse of the bubble, which results in a 
small difference at the base of the jet where the free surface connects to the wall of the 
capillary. This however has no significant effect on the part of the jet that is further away 
from the contact point. 

3.1. Jet and bubble shape 

Figure [4] shows how the jet and the bubble develop in time. The maximum speed of the 
jet is approximately 30 m/s, which is reached around 20 /is after the pressure pulse. The 
jet has a diameter which is about 1/10 of the diameter of the tube; this holds for all tube 
diameters that we tested. The bubble initially grows spherically but, due to confinement 
and asymmetry, it later takes on an elongated shape and grows more towards the free 
surface. The right side of the bubble is almost stationary, also during the collapse where 
a thick jet is formed reminiscent of the collapse of a bubble near a free surface or a solid 
boundary (Blake & Gibsonjl981; Lindau fc Lauterborn|2003 |. 



3.2. Jet velocity and velocity field 

In this paper we focus on the jet velocity which can directly be compared with the 
experimental measurements. We define the jet velocity as the velocity component parallel 
to the tube axis at the tip of the jet (in figure [4] to the left). Figure [5] shows numerical 
results together with experimental results on how the jet velocity evolves in time. There 
are two acceleration mechanisms: First, driven by the very short pressure pulse, a speed 
of about 13 m/s is reached almost instantaneously. After this there is no more driving, 
but the focusing of the flow accounts for a further acceleration of the jet which reaches 
a maximum velocity of about 25 m/s. Deceleration is accounted for by surface tension 
(see 



3.4 1 and the collapsing bubble. 



Figure 



6] shows the velocity field in the liquid during jet formation. In figure [6fa) the 
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Figure 4. The evolution of the jet and the bubble from a boundary integral simulation. 
Parameters: Rtube = 250 /im, Ap = 2027 bar, At = 50 ns, 6 — 30 degrees, A = 1106 /im. 
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Figure 5. The velocity of the tip of the jet as a function of time for both the simulation (solid 
line) and experiments (diamonds). After an almost instantaneous acceleration to 13 m/s during 
the 50 ns pressure pulse, the tip is further accelerated by the focusing geometry to about 25 m/s. 
Experimental conditions and numerical settings are the same as in Fig. [I] 



interface has not moved significantly due to the small time interval, but it clearly shows 
how the interface is responsible for the focusing of the flow. In figure [6F&) we see that 
although the surface in the center is approximately flat, the velocity still has a focusing 
profile. Indeed, at t = 7.5 ^<s the jet is still accelerating. Only after about 15 /is (see 
figure |6[e)), there is no more focusing of the flow. Stretching of the jet is visible in 
figure |6[rf) and[6]je) where the velocity of the tip of the jet is larger than the velocity at 
the base, with the consequence that the jet is becomes thinner while increasing in length 
( |Eggers fc Villermaux|[2fJ08l ). 



3.3. Contact angle 

Now, how does the jet speed depend on the contact angle? As we explained in 
meniscus initially has the shape of a spherical cap. The meniscus thus has a well-defined 
radius of curvature, which depends on the contact angle 8. A contact angle of 90° results 



2.1 the 
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Figure 6. The velocity field during the formation of the jet. (a) is right after the pressure pulse, 
the free surface still has a spherical cap shape. In (a), (b) and (c), focusing of the flow can be 
seen, (d) and (e) clearly show the stretching of the jet: The largest velocity is in the tip of the 
jet, and gradually decreases towards the base of the jet. 
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Figure 7. The influence of the contact angle on the jet velocity. Tube radius is 250 fxxa, and 
A = 460 pm. The numerical results are obtained with Ap = 1647 bar, the absorb ed laser energy 
in the experimental data was 458 /xj. The dashed line represents the model ( |4.14[ ), with a = 0.44, 
/3 = 1.33, and ho = 0.26. 



in a fiat interface (no curvature), and a contact angle of 0° gives a radius of curvature 
equal to the inner radius of the capillary. A smaller contact angle increases the curvature 
of the free surface, and therefore increases the focusing of the flow. 

In the experiments, 9 can be measured directly from images of the static meniscus, so 
we can directly compare the influence of the contact angle in experiments and simulations. 
Figure [7] shows that the jet velocity clearly increases with cos 9. A more careful look 
reveals that the dependence on the contact angle is a bit stronger than linear, both in 
the experiment and the computation. There is a non-linear dependence on cos 9 because 
the curved interface increases the velocity in two stadia of the jet formation: During the 
pressure pulse and during the flow focusing. In §H] we will provide a detailed analysis 
of this dependence on 9. The absorbed laser energy in the experiments in figure [7J was 
458 /zJ, the applied pressure in the simulations was 1674 bar. 



3.4. Surface tension 

In order to study the effect of surface tension in isolation, we turn off the collapse mech- 
anism of the bubble by setting the pressure in it to atmospheric after the initial pressure 
pulse. In this case the bubble keeps growing until it would be ultimately restrained by 
surface tension over a much longer time scale than we consider. 

Figure [8] shows the development of the jet velocity for four different values of the 
surface tension. Clearly, the acceleration phase is dominated by inertia, as there is almost 
no difference in the acceleration while there is an order of magnitude difference in the 
surface tension. Only when the jet reaches its maximum velocity and during deceleration 
surface tension starts to play a role. This stands to reason, because the only decelerating 
mechanism in this case is surface tension, counteracting the increase of surface area 
caused by the jet. 

We can calculate the order of magnitude of the expected decrease in velocity V s using 
an argument similar to the Taylor-Culick velocity of a liquid sheet ( Taylor(l959 Culick 
1960 Gordillo & Gekle 2010): Using the jet radius Rj as the relevant length scale, we 
can balance the kinetic energy per unit length E k ~ ^npR'jV^ with the surface energy 
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Figure 8. Jet speed versus time for different surface tensions. The solid line has the surface 
tension of water (a = 72.8 mN/m), the other values for the surface tension are a = 35 mN/m 
(dash-dotted line), a — 150 mN/m (dashed line) and a = 300 mN/m (dotted line). The collapse 
of the bubble was turned off in these simulations to isolate the effect of surface tension. Tube 
radius is 250 fim, A = 1106 /im, and Ap = 2027 bar. For reference, the evolution of the jet 
velocity with bubble collapse is represented by the gray solid line, with the same parameters 
and a = 72.8 mN/m. 



of the jet per unit length E s ~ 2iraRj and arrive at a velocity 



Vs 



4(7 

pRj 



1/2 



(3-1) 



Using Rj m 25 /um and the material properties of water gives V s ~ 3.4 m/s, which is 
indeed similar to the decrease in velocity presented by the solid line in Figure [8] 

3.5. Pressure impulse 

The magnitude of the pressure impulse ApAt is the only parameter that cannot be 
directly related to the experiment. We expect that the absorbed laser energy is the 
experimental parameter that is most directly related to it. In the experiments, a linear 
relation between the energy and the velocity of the jet is found. As can be seen in figure[9j 
in the numerical simulation the jet velocity also depends linearly on the pressure pulse. 
From these observations we conclude that there is a linear relation between the absorbed 
laser energy (E in Joules) as measured in the experiments and the pressure impulse that 
we apply in the simulations of figure [9] 

ApAt « 3.30- 10 4 sm" 3 £- 5.0 Pa- s (3.2) 



for Rtube — 250 pm, and 



ApAt sw 4.39 • 10 5 s m~ 3 £-8.5Pa-s 



(3.3) 



for Rtube — 25 /im. The prefactor is one order of magnitude larger when the tube radius 
is one order smaller, while the threshold value is of the same order. More about the 
influence of the tube radius can be found in § |3.6| 

In both the simulation and the experiment, there is an apparent threshold value for the 
energy or pressure below which we can not observe a well-defined jet. In the experiments, 
the main reason for this would be that a large amount of the laser energy is lost in heating 



up the fluid before a bubble can be created, as was shown by Sun et al. (20091. In the 
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Figure 9. The influence of the pressure on the jet velocity, compared to experiments for two 
different tube radii, (a): Tube radius 250 /im, At = 50 ns, 6 — 30 degrees and A = 600 /im. 
(6): Tube radius 25 /im, At = 5 ns, 9 — 30 degrees and A = 400 /im. Experimental data were 
converted from the absorbed laser energy to pressure impulse by a fitting routine; the energy is 
indicated at the top axis. 

simulations the only cause for the threshold lies in the surface tension that prevents the 
formation of a jet, and the vacuum inside the bubble after the pressure pulse is applied. 
Once the kinetic energy is much larger than the surface energy related to the formation 
of the jet and the potential energy related to the size of the vacuum bubble, a jet can be 
formed. The zero value of E extrapolated on the upper horizontal scale in figure [9] lies 
considerably to the left of the zero value of the lower impulse scale, which implies that 
the experimental threshold, due to thermal and other losses, is significantly higher than 
the numerical one. 



3.6. Tube radius 

In the experiments, there is a clear dependence of the jet velocity on the tube radius 
Rt- smaller tubes create faster jets with the same absorbed laser energy, approximately 
following the relation Vj et oc 1/Rt- One naively would argue that a smaller tube will 
provide a stronger curved free surface, and therefore the stronger focusing will result in a 



faster jet. In contrast, Figure 10 shows that the maximum jet speed that is obtained in the 
simulations at most only shows a very weak dependence on the tube radius (if there is any 
at all). However, a different effect caused by the tube radius is very clear: The acceleration 
is much larger for smaller tubes, so that the maximum velocity is reached earlier. Indeed, a 
smaller tube provides a higher curvature of the free surface, so the acceleration due to flow 
focusing is larger. The maximum velocity, however, is a combination of the magnitude and 
the duration of the acceleration, which both depend on the tube radius. The simulations 
show that these two parameters cancel each other if we only change the tube radius: larger 
tubes have less acceleration due to focusing of the flow, but the acceleration persists for 
a longer time, as can be seen clearly in figure [TUj This is consistent with dimensional 
analysis, provided that the relevant time and velocity scales are taken to be Rt/Vo and 
Vq = ApAt/(pX), which will be motivated and discussed in §E] This scaling will make 
the curves collapse approximately up to the maximum, confirming that the tube radius 
does not influence the jet velocity. After the maximum there is no collapse due to the 
influence of surface tension (cf. § 3.4). 
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Figure 10. The influence of the tube radius on the jet velocity. A smaller radius increases the 
acceleration due to flow focusing, but decreases the duration of the acceleration, resulting in 
approximately the same maximum jet velocity. Tube radii in this figure are 25 /tm (solid line), 
50 fim (dashed line) and 100 fim (dotted line). The other parameters are the same for all three 
simulations: Ap = 3040 bar, At = 50 ns, 9 = 30 degrees and A = 1250 /^m. 



Although our simulations and analysis show that the jet velocity is independent of the 
tube radius, the remaining question is why there is such a strong dependence on the tube 
radius in the experiments. The most plausible explanation is that for a fixed absorbed 
energy, the generated pressure has a strong dependence on the tube radius. The reason 
for this could be found in the volume V e that is exposed to the laser energy E. Based on 
dimensional analysis the produced pressure can be expected to scale as 



Ap ex 



E 

Ve. 



(3.4) 



A smaller tube would result in a smaller volume that is exposed to the laser energy. This 
then would account for the dependence of the jet velocity on the tube radius that was 
observed in the experiments. 



3.7. Bubble distance from the free surface 

In experiments, the most direct measurable parameter is the distance between the menis- 
cus and the bubble. This makes it an excellent candidate to compare with numerical 
simulations. We define the distance A as the distance between the center of the bubble 
and the point on the meniscus that is on the axis of symmetry (see figure [2]) . Due to 
the axisymmetry of the numerical simulations, the bubble is always in the center of the 
tube. In experiments however, the bubble is usually created near the wall of the tube, 
due to the characteristics of the absorption of the laser light in the liquid. The difference 
between these different bubble positions can be neglected when the distance between the 
bubble and the meniscus is large compared to the radius of the tube (X/R t ~> 1). 

Figure [TTJ shows the jet velocity as a function of A, together with the experimental 
measurements. There is a good agreement between experiments and simulations, and 
both show a clear 1/A dependence for the jet velocity. Both the numerical and the 
experimental results are obtained with a tube radius of 250 fim and a contact angle of 
25 degrees. The absorbed laser energy in the experiments was 305 /LtJ and the pressure 



amplitude for the simulations, calculated using (3.2 1, was 1013 bars 
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Figure 11. The influence of the distance between the bubble and the meniscus on the jet veloc- 
ity. Tube radius is 250 /im and 6 = 25°. The numerical results are obtained with Ap = 1013 bar, 
directly related with equation ( |3.2| ) to the absorbed laser energy in the experimental data, which 
was 305 /iJ. A power-law fit (solid line) reveals a clear 1/A dependence of the jet velocity. 



4. Theoretical approximation 

We will now try to analytically understand the mechanism of the jet formation and 
the achieved velocities by approaching the problem with a simplified model. The route 
through which the jet attains its maximum velocity can be split in two parts: the pressure 
impulse and the flow focusing. The effect of the former is determined by how the pressure 
wave, which gives an initial velocity to the liquid, is modelled in the (incompressible) 
simulations which will be discussed in § |4.1| After the liquid is set into motion, the 
curved shape of the meniscus leads to a further acceleration of the liquid by focusing it 



in a fast thin jet (§ 4.2 1. We first neglect the influence of the curvature of the free surface 



on the action of the pressure pulse to compute the maximum jet velocity (§4.3), and 
then revisit this approximation in §|4.4| 



(4-1) 



4.1. The pressure pulse 

We approximate the system during the pressure pulse as one-dimensional, so after ne- 
glecting viscosity and compressibility we can write the Euler equation as 

du 1 dp 

dt pdz 7 

where u is the axial velocity of the liquid, and t the time. Note that due to continuity in 
this one-dimensional system, the du/dz term in the material derivative equals 0, so on 
the left hand side we only have du/dt. The axial pressure gradient dp/dz is given by the 
pressure in the bubble Ap and the distance A between the bubble and the free surface: 

dp Ap 

dz = X' 

The Laplace pressure jump on the free surface can be neglected because Ap is very large 
compared to the typical pressure associated with surface tension. 

We integrate |4 . 1 1 over the duration At of the pressure pulse and obtain the velocity Vo 
of the free surface after the pressure pulse ( Ory et alj|2000 ): 

ApAt 



(4.2) 



Vq = 



(4.3) 
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where we assume A to be constant, which can be done if At is small enough. 

4.2. The flow focusing 

After the pressure pulse there is no more driving of the flow, which means that all further 
acceleration is caused focusing. 

We will now give an analysis for the acceleration due to flow focusing based on con- 
tinuity. Starting with a spherical surface with radius of curvature R c and velocity Vq 
directed normal it, we keep the flow rate constant: 

V R 2 C = (V + dV)(R c - dR) 2 , (4.4) 

where dV is a small increase in velocity due to a small decrease in radius dR. At leading 



order, dR = V dt, and (4.4 1 becomes 



f - ¥ 

The radius of curvature can be expressed using the tube radius Rt and the contact angle 
8 as R c = R t / cos9, which then gives us the following expression for the acceleration: 

a = 2Vt^ (4.6) 

Note that the same scaling for the acceleration can be obtained by dimensional analysis 
using Vq and R c as the relevant velocity and length scale respectively. 

Clearly, smaller tubes have stronger focusing and therefore generate a larger accelera- 
tion. This, however, does not mean that the maximum velocity of the jet will be higher 
as well. To determine the increase in speed due to the flow focusing we have to find a 
time scale during which the fluid is accelerated. The focusing time scale Atf is provided 
by the typical velocity Vo (the velocity created by the pressure pulse) and the typical 
length scale R t (the radius of the tube) : 

At f = ^. (4.7) 

The increase in velocity due to flow focusing is then: 

AV ~ aAt f = 2V Q cos 9, (4.8) 

where it becomes clear that the increase in velocity due to flow focusing is independent 
of the tube radius. 

4.3. The maximum jet velocity 

The maximum velocity reached by the jet is the sum of the velocity reached after the 
pressure pulse and the increase in velocity due to flow focusing: 

V max = V Q + AV=^±(l + [5cos6) (4.9) 

p\ 

with /3 a proportionality factor which we expect to be of order unity. First of all, the 
proportionality to p x l is in excellent agreement with the results from the simulations 
shown in figures M fij] and with the fact that V max does not depend on the tube radius 



(see figure 10). To compare the dependence on the curvature of the meniscus, we now 
turn to figure [T3|( a). Here we compare the model with (3 = 2.0 to the simulation data. 
The velocities are roughly reproduced, but it is clear that there is a dependence on the 
curvature for Vq, which is not accounted for by the model and the increase of V max is 



Highly focused supersonic microjets 



15 




Figure 12. Definition of the distance Ho, the effective distance H e 

and angle f. 



radius of curvature R c 



therefore not very accurately reproduced. Clearly, neglecting the curvature of the surface 
during the pressure pulse has been too bold an assumption. 

4.4. Correction for a pressure pulse on a curved interface 
We will now apply a correction to the above derived model to account for the curved 



interface during the pressure pulse. The above derivation (4.3) gives the velocity in the 



bulk, far away from the the bubble and the free surface. We will use volume conservation 
and an approximate velocity distribution on the free surface to calculate the free surface 
velocity on the tube axis. 

The first step will be to determine the velocity distribution on the free surface. Due 
to the short time scale and the magnitude of the pressure pulse, we can neglect the 
tangential velocity components. We therefore only take into account the velocity normal 
to the interface, and we can consider the free surface as an equipotential surface. Away 
from the free surface we expect the one-dimensional approximation to hold, resulting in a 
uniform axial velocity, so that we will have evenly spaced equipotential surfaces oriented 
perpendicular to the tube wall. 

The free surface is a curved equipotential surface that needs to matched to the plane 



equipotential surfaces in the bulk. With reference to figure [12] we introduce a distance 
Hq, ultimately to be treated as a fitting parameter, as the smallest distance from the 
free surface where we assume the equipotential surfaces to be unaffected by the curved 
interface. We now calculate the distance between the free surface and the horizontal 
plane defined by H in the direction normal to the free surface. We will call this the 
(■-dependent effective distance H e : 



He(0 = 



R c + Ho 
cos( 



-R« 



(4.10) 



where we have defined the position on the free surface as a function of the radius of 
curvature R c and the angle (, as shown in figure |12| Because the potential difference 
between the plane defined by Hq and the free surface is constant, we expect the velocity 
to be inversely proportional to the effective distance H e : 

VMO = ^ (4.ii) 

with A a constant that we will determine using volume conservation: The flux through a 
cross-section in the bulk, where the velocity is uniform, must equal the flux through the 
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Figure 13. The velocity Vb of the jet just after the pressure pulse and the maximum velocity 
Vmax of the jet. Parameters: R tu be = 100 /im, Ap = 3040 bar, At = 50 ns, A = 1250 /mi. 
The diamonds and circl es a re results from boundary-integral simulatio ns. In (a) the solid line 
corresponds to equation 4.3 and the dashed line corresponds to equation 4.9 with j3 = 2. In (b) 
the solid line corresponds to equation |4.13| and the dashed line corresponds to equation |4.14[ 
with a = 0.94, /3 = 0.84, and h = 0.26. 



free surface: 



V B TTRf = 



2ir /•7T/2- 



V fs (()R 2 c sin (d(d(;, 



(4.12) 



where Vb 



>o Jo 
i p x l , R c = R t j cos 6, and R t the tube radius, a is aprefactor which should 
be of order unity, reflecting the one-dimensional character of the flow in the bulk. We 
expect a to become closer to 1 when X/Rt increases. 



We now have an expression for A, which we substitute in (|4.11|), and we arrive at the 
following velocity on the free surface Vq 



Vfstt 



V 



ApAt 1 



= 0,t = 

cos 6 



At): 



pX 2h obl 



( sin9-b \ 



■sin6>- 1 



(4.13) 



with b = 1 + h cos 8, and the geometrical factor h a = H a /R t . The value of ho only 
needs to be determined once by fitting, since we do not expect it to change with other 
parameters. 

The maximum velocity remains 



V maB = V (l + pcoa9), 



(4.14) 



with Vb given by (4.13). 



In figure |l3[ fc) we compare the corrected model to the boundary integral simulations, 
with a = 0.94, (3 = 0.88, and ho — 0.26. We find an excellent agreement between the 
simulations and the model given by (4.131 and (4.14). Figure [7] shows the comparison of 
the model with both experiments and simulations, for a different tube radius and bubble 
distance, but we have used the same value for ho- Note that in the comparison shown in 
figure [7J a is smaller, which is due to the fact that in that case X^ Rt- 



5. Conclusions and discussion 

We have numerically investigated the formation of microjets in a capillary by laser 
induced cavitation using axisymmetric boundary integral simulations. Although com- 
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pressibility plays an important role in the formation and initial growth of the bubble as 
well as in the subsequent pressure wave that travels through the liquid, we have assumed 
incomprcssibility for our numerical simulations. This is possible because the compress- 
ibility is only important during the very short period of the pressure wave reflecting from 
the free surface, which we have modeled by applying a short strong pressure pulse on 
the bubble inside the capillary. After the initial pressure impulse, the formation of the 
jet can be considered as incompressible, because the observed speeds are mostly much 
smaller than the speed of sound in water, and pressures are moderate. 

We have found a convincing agreement in shape and evolution of the jet between the 
simulations and the experiments, which has allowed us to perform a detailed study of 
the involved parameters, including those which are difficult to access in experiments. 

We compared the influence of the different parameters on the maximum achieved 
velocity of the jet. We have found good agreement between the simulations and the 
experiments by investigating the influence of the distance A and contact angle 9. It is 
however much harder to directly relate the absorbed laser energy in the experiment to the 
pressure pulse in the simulation. By comparison we were able to derive that the pressure 
pulse ApAt is linearly related to the absorbed energy E, and we have given arguments 
for the dependence of the pressure on the capillary radius R t . 

The jet velocities we find in the simulations can be reproduced accurately by a simple 
model. We developed this model starting with a one-dimensional approximation for the 
pressure pulse and dimensional analysis for the focusing effect during jet formation. 
We improved the one-dimensional approximation by making a correction for the curved 
interface during the short pressure pulse, where deformation of the meniscus can be 
neglected. 

We would like to thank Claas Willem Visser, Sascha Hilgenfeldt and Leen van Wijn- 
gaarden for helpful discussions. LP. especially thanks Stephan Gekle for his assistance 
with the BI code, and acknowledges NWO for financial support. 
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